Spin Accumulation in Quantum Wires with Strong Rashba Spin-Orbit Coupling 
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We present analytical and numerical results for the effect of Rashba spin-orbit coupling on band 
structure, transport, and interaction effects in quantum wires when the spin precession length is 
comparable to the wire width. The situation with only the lowest spin-split subbands occupied 
is particularly interesting because electrons close to Fermi points of the same chirality can have 
approximately parallel spins. We discuss consequences for spin-dependent transport and effective 
Tomonaga-Luttinger descriptions of interactions in the quantum wire. 
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Spin-dependent transport phenomena are currently at- 
tracting a lot of interest because of their potential for fu- 
ture electronic device applications 0. Basic design pro- 
posals for spin-controlled field-effect switches!^ ^ use the 
fact that electron waves with opposite spin aquire differ- 
ent phase factors during their propagation in the pres- 
ence of Rashba spin-orbit coupling^ (RSOC). The latter 
arises due to structural inversion asymmetry in quantum 
heterostructures||, |^ where two-dimensional (2D) elec- 
tron systems are realized. The single-electron Hamilto- 
nian is then of the formQ] H21) = Hq + Hso where 
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with m denoting the effective electron mass[g^. Possibil- 
ity to tune the strength of RSOC, measured here in terms 
of the characteristic wave vector fcgo, by external gate 
voltages has been demonstrated experimentally!^, ||, p^ . 
As a manifestation of broken spin-rotational invariance, 
eigenstates of i?2D which are labeled by a 2D wave vec- 
tor k have their spin pointing in the direction perpen- 
dicular to k. Hence, no common spin quantization axis 
can be defined for eigenstates when spin-orbit coupling 
is present. Confining the 2D electrons further to form a 
quantum wire, one might naively expect to again be able 
to define a global spin quantization axis, as the propa- 
gation direction of electrons in a one-dimensional (ID) 
system is fixed. However, this turns out to be correct 
only for a truly ID electron system with vanishing width. 
In real quantum wires, such a situation is approximately 
realized when the spin- precession lengthj^ Tr/fcgo is much 
larger than wire width. Another way to formulate this 
condition is to say that the characteristic energy scale 
Ago = fi-^^so/^™ for RSOC is small compared to the en- 
ergy spacing of ID subbands. For a quantum wire defined 
by a parabolic confining potential, e.g., 
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the latter would be Tiuj. When spin-orbit coupling is 
not small (i.e., when Ago ~ fii-^ for the case of parabolic 



confinement), hybridization of ID subbands for opposite 
spins becomes important, resulting in the deformation of 
electronic dispersion relations jl^ . The effect of this de- 
formation on transport properties has been the subject 
of recent investigation p]| , e.g., with respect to implica- 
tions for the modulation of spin-polarized conductances 
as a function of RSOC strengthjl^ which is the principle 
of operation for spin-controlled field-effect devices]^, ^. 

Here we present results for the detailed spin struc- 
ture of electron states in a quantum wire, defined by 
the parabolic confining potential V{x) given in Eq. (|^), 
with strong RSOC present. Contrary to previous as- 
sumptions that were uncritically adopted in the recent 
literature |lj], we find that electrons with large wave vec- 
tors in the lowest spin-split subbands have essentially 
parallel spin. The spin state that right-moving electrons 
converge toward is opposite to that for left-movers. This 
counterintuitive result will be explained qualitatively in 
the following paragraph, before presenting analytical and 
numerical results for electronic dispersion curves and spin 
structure of eigenstates. A texture-like variation of spin 
density across the wire is identified. We then apply 
the Landauer-Biittiker formalism |l^, |l^ to discuss spin- 
dependent transport in hybrid systems of a wire with 
RSOC attached to leads where kso — 0. Current turns 
out to be spin-polarized in the wire but unpolarized in 
the leads. We elucidate the peculiar current conversion 
at wire-lead interfaces that sustains this novel type of 
spin accumulation in the wire. Finally, consequences for 
the low-energy description of interacting wires in terms 
of Tomonaga-Luttinger-type models are discussed. 

We start by considering basic features for eigenstates 
of the Hamiltonian _ffiD = i?2D + which are ID 

plane waves in the y coordinate with wave number ky 
but bound in x direction. At finite kgo, spin degeneracy is 
preserved only for eigenstates with ky = 0; their energies 

are the shifted harmonic-oscillator levels i^i"' — ^(2n + 
1) — Ago. This result is exact. To characterize states with 
finite ky , we rewrite i?iD = ^^pb + ^^mix where 
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FIG. 1: Lowest and first excited spin-split subbands of a quan- 
tum wire, defined by a parabolic confining potential with os- 
cillator length in a 2D electron system, with strong Rashba 
spin-orbit coupling such that kso luj = 0.9. Thick curves are 
results of the exact numerical calculation, while thin curves 
are obtained using the approximate two-band model which 
includes only spin-orbit induced mixing of the lowest two 
parabolic subbands. Evidently, this approximation gives rea- 
sonable results for the lowest spin-split subband, even in the 
present case of a rather large spin-orbit coupling strength. 
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and i/mix — —TT-ksoOyPx/m. Straightforward calculation 
yields eigenstates of i?pb which are also eigenstates of (Jx 
with eigenvalue cr = ±1 and have energies E^J^\ky) = 
%^(2n + 1) + |i {ky + (Tfcso)' - Aso. The term H^^^ 
induces mixing between the shifted parabolic subbands 
Ei^J^\ky). To lowest order in perturbation theory, it re- 
sults in a uniform shift of eigenenergies by —Ago and a 
small deviation of spin quantization in x direction [pT|. 
Hence, for Ago ^ eigenstates of 7?id have ener- 
gies E^a\ky) — Ago and are, to a good approximation, 
eigenstates of ax ■ When Ago becomes comparable to the 
subband splitting, anticrossings occur between neighbor- 
ing subbands with opposite spin index a. As a result, 
no common spin-quantization axis can be defined any- 
more for eigenstates within any subband. Far enough 
from anticrossings, eigenstates of Hso will essentially be 
eigenstates of ffpb- In particular, their spins will be ap- 
proximately aligned in x direction. In the lowest two sub- 
bands, right-movers with wave vectors larger than that of 
the anticrossing point can then have approximately par- 
allel spin. The same is true for left-movers whose asymp- 
totic spin direction is opposite to that of right-movers. 

In Fig. |lj we show as thick lines numerically calcu- 
lated spectra of Hyd for a large value of spin-orbit cou- 
pling. Deviation from parabolicity is clearly visible. In- 
terestingly, it is possible to obtain a good quantitative 
description of the lowest spin-split subband by diagonal- 
izing Hyd in a truncated Hilbert space which is spanned 
by the lowest and first-excited spin-degenerate parabolic 
subbands of the Hamiltonian H[) + V{x). We call this the 
two-hand model and find an approximate expression for 



FIG. 2: Spin structure of electron states in a quantum wire 
with strong spin-orbit coupling. Panel a): Expectation value 
of spin projection onto the x direction for electron states ob- 
tained in Fig. |l|. Results of exact numerical calculation for the 
lowest spin-split subbands (main figure) and first excited spin- 
split subbands (inset) are given by thick curves. The effective 
two-band model reasonably approximates the behavior of the 
lowest subband (thin lines in the main figure). Right-moving 
electrons with large wave vectors asymptotically have parallel 
spin which is opposite to that of left-movers. The same can be 
observed in panel b) where the spectrum in a finite magnetic 
field B pointing in x direction is compared with that in zero 
field. Here dispersion curves are calculated within the two- 
band model, for Zeeman energy gn^B = Q.25huj (thin lines) 
and in zero magnetic field (thick lines). 

the dispersion of the lowest spin-split subband, 
2^(2b) 

= 2 + {kyl^f - - 72fcsofc,/2 )2 + 2(fc,oL)2 , 

(4) 

where = y/fijmuj is the oscillator length of the 
parabolic confinement, and 7 = ± a subband index that 
does not have the meaning of a spin quantization num- 
ber. We show Eq. (Q) and the corresponding result for 
the first excited subband as thin lines in Fig. |l|. It is seen 
that the two-band model is quite adequate for the lowest 
subbands, even for rather strong spin-orbit coupling. 

Results shown in Fig. ^ confirm conclusions reached 
in our previous discussion of the spin structure of elec- 
tron eigenstates with RSOC present. Panel a) shows the 
expectation value of spin component in x direction for 
eigenstates of -ffiD in the lowest and first excited spin- 
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FIG. 3: Texture- like structure of spin density across the quan- 
tum wire, calculated within the two-band model for states in- 
dicated by black dots in Fig. ^, which have energy 0.75?J.tJ. 
a) Spatial variation of nonzero components of spin density for 
the state with larger wave vector, b) Same for the other state, 
c) Visualization of spin texture for the same state as a). Ar- 
row length is proportional to spin density, d) Spin texture 
visualized for the same state as considered in b). 



split subbands for the same value of fcgo used in Fig. |^. 
Data in Figs. ^ and || for the same subband are indi- 
cated by the same line type. For the lowest subbands, 
we also give, as thin lines, results obtained analytically 
within the two-band model. It is clearly seen that spins 
of eigenstates with large absolute value of wave num- 
ber are approximately quantized in x direction, with spin 
direction of left-movers being opposite to that of right- 
movers . This fact is underscored by the properties of 
the energy spectrum in a finite magnetic field B va. x di- 
rection which is shown in panel b). Clearly, Zeeman shift 
of states at large positive wave number is opposite to 
that for states with large negative wave number. Shown 
as thin lines in the main figure of panel a) are curves 
obtained analytically within the two-band model which 
yields again reliable results for the lowest subbands. We 
therefore use it to calculate the variation of spin density 
s[x) = '^'^ {x)d'^(x) across the wire. [The spinor ^{x) 
denotes the transverse part of an eigenfunction of i?iD 
which, in the presence of spin-orbit coupling, depends 
on wave vector.] It turns out that the density Sy{x) of 
spin component parallel to the wire vanishes identically. 
Hence, only the x and z components of spin density are 
shown in Fig. |^, displaying an interesting texture-like 
variation with coordinate x whose structure reflects the 
mixing between subbands due to -ffmix- Note that the 
expectation value for the z component of spin vanishes 
for eigenstates of i?iD- 



From the above it has become clear that, in general, 
spin quantum number is not an appropriate way to char- 
acterize electron states in a quantum wire with strong 
RSOC. Only states with wave number ky far enough from 
anticrossing points will asymptotically have their spin 
quantized in x direction. From considering Figs. |l| and 
^, the following special situation can be envisioned which 
has rather counterintuitive consequences. At low enough 
electron density such that only states in the lowest spin- 
split subbands are occupied, states near the Fermi en- 
ergy sf will be localized near four Fermi points. When 
electron density is not too low, their spins are approxi- 
mately quantized in x direction. As pointed out above, 
spins of states near Fermi points for right-movers are ap- 
proximately spin-down, opposite to the spin direction of 
left-moving states near ep. Assuming it to be possible 
to selectively raise (lower) the electrochemical potential 
of right-movers (left-movers), a spin-polarized current 
could be generated. Usually, creating a population of 
left-movers and right-movers with different electrochem- 
ical potentials is achieved by coupling the quantum wire 
adiabatically to ideal contacts||l^, However, the un- 
derlying assumption that excess electrons injected from 
the right (left) reservoir will only be spin-up (spin-down) 
is not realistic because, typically, RSOC will be absent in 
the contacts. The different nature of electron states in the 
wire and the leads will result in strong scattering at wire- 
lead interfaces. Similar to the approach taken in Ref. [ p2[ , 
we model this situation by attaching semi-infinite leads 
with kso = to the wire where fcgo 7^ 0. The transmission 
problem can be solved exactly by matching appropriate 
ansdtze for wave functions in the wire and the leads. The 
usual condition for ensuring current conservation has to 
be modified because the group velocity for electrons in 
the quantum wire with RSOC readsjlj, |l^ 

Vy = h{ky + kso(Tx)/m . (5) 

Despite the unusual spin structure at the four Fermi 
points which is asymmetric with respect to right-movers 
and left-movers, no spin-polarized current is generated 
in the leads. However, as is shown in Fig. ^ a process 
of current conversion occurs close to the interfaces in the 
wire that results in a finite spin polarization of current in 
the wire. We have therefore found a unique type of spin 
accumulation that is not, as in the usual case|2^ , induced 
by ferromagnetic contacts. Our analysis shows that cur- 
rent conversion is enabled by scattering into evanescent 
modes of the wire because of the peculiar form of the ve- 
locity operator (^. A four-terminal measurement with 
ferromagnetic contacts as weakly coupled voltage probes 
should enable experimental verification of spin accumu- 
lation in the wire. 

Finally, we briefly remark on the effective low- 
energy description of an interacting quantum wire with 
strong RSOC. In the spirit of Tomonaga-Luttinger 



models pil, 22 for interacting ID systems, we linearize 



4 




5 10 







,' 1 / 

1 1 


1 ' 


r '\ 


1 

.■\ '\ 




! 1 1 
I ' ' 


\, w 


WW 

li 

It 


' 1 ( I , 

" ' 


I 

I / 1 
I / ' 
\/ I 

1 




1 


1 


1 





5 10 15 20 



FIG. 4: Transport in hybrid systems of a wire with strong 
Rashba spin-orbit coupling and ideal leads, calculated exactly 
using the Landauer-Biittiker formalism within the two-band 
model (inclusion of higher subbands leads only to small quan- 
titative changes). Panel a) shows the spatial variation of cur- 
rent polarization in a semi-infinite wire {y > 0) attached to 
an ideal lead {y < 0). Conversion of incident spin-up current 
is illustrated in the inset. Here 7o-| denotes the spin-a current 
in the wire when spin-| current is injected from the lead. A 
finite spin polarization exists also in a finite wire with to semi- 
infinite leads attached [panel b)]. Here quantum interference 
gives rise to additional oscillatory structure. Parameters used 
in the calculation are Ep = 1 hui and kso = 0.9 / J^ 

the single-electron energy spectrum close to the four 
Fermi points. We expUcitly avoid attaching any spin la- 
bels. Rather, we define type-A (type B) right-movers and 
left-movers having the same velocity wa (''^b)- Typical 
electron-electron interactions give rise to a term Hint — 
i / x.y 'ip^tp{x,y)U{x — x',y — y')'tp'''ip{x' ,y') in the elec- 

x' .y' 

tron Hamiltonian. In the low-energy, long-wave-length 
limit, we can write ii{x,y) = ^ „=a,b V'a/3(2/)^'fcF„f) (2;) 
and assume U to be long-range on the scale of the wire 
width but short-range on the scale of the wire length. 
It is important to note that the present case differs from 
the usual one in that the transverse wave- function spinors 
(x) are nearly orthogonal. As a result, backscatter- 
ing processes are strongly suppressed. Apart from this 
fact and the peculiar spin structure of states near the 
four Fermi points, the present system is identical, on a 
formal level, to a two-component |p3[ or Zeeman-split|24 
Tomonaga-Luttinger model. Response to an external 
magnetic field will, however, be special in the present 
case. Postponing a detailed analysis to a later publi- 
cation, we mention here only a few basic facts. When 



Fermi points are far enough away from anticrossings, a 
magnetic field B applied in x direction will shift right- 
movers (left-movers) to higher (lower) energies. [See 
Fig. ^).] The Zeeman term in bosonized form reads 
then Hz — Ix^P' "^tiere Hp is canonically con- 

jugate to the phase field 6p that is related, within the 
usual [^ phase-field formalism, to the total electron den- 
sity ptot = I]°=^.B Pq/3 via y^2jTrdy0{y) = ptot(2/)- Ap- 
proximate orthogonality of transverse parts of electron 
wave functions enables spin-flip processes, in the long- 
wave-length limit, only between left-moving and right- 
moving branches of the same type (A or B). In general, 
any spin-flip process incurs a large momentum transfer. 
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